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ABSTRACT 

By assuming a particular mass function we find new exact solutions to the Einstein 
field equations with an anisotropic matter distribution. The solutions are shown to be 
relevant for the description of compact stars. A distinguishing feature of this class of 
solutions is that they admit a linear equation of state which can be applied to strange 
stars with quark matter. 
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1 INTRODUCTION 

In studies of compact relativistic astrophysical objects, 
strange stars have long been proposed as an alternative to 
neutron star models. The theoretical basis of strange stars 
composed of deconfined u, d and s quarks is a direct con- 
sequence of the conjecture that quark matte r might be the 
true ground state of hadrons (|WitterJ[T98i : iFarhi fc Jaffel 
1984). Strange stars are expected to form during the col- 
lapse of the core of a massive star after the supernova explo- 
sion. Another possibility is that a rapidly spinning neutron 
star can accrete sufficient mass to undergo a phase tran- 
sition to become a strange star. Although we have no di- 
rect experimental evidence which may support the existence 
of such objects, recent observational data on the compact- 
ness (mass-radius ratio) of some of the compact objects like 
Her X-l, SAX J1808.4-3658, RX J185635-3754, 4U 1820-30, 
4U 1728-34, PSR 0943+10, strongly favour the possibility 
that they could act ua lly be strange s tars | Bombaci| 1 19971 ; 
Li et aJll995l. Il999al lbl: iDev et alliwi ; IXu et aj|l999l l200ll ; 



Pons et a i 2002). Naturally, the internal composition and 



consequent geometry of such objects have become a subject 
of considerable interest over the last couple of decades. As 
the physics of very high density matter is still not very clear, 
most of the strange star studies have bee n performed within 
the fr a mework of a bag mo del, see e.g., iGlendenning et ~a\ 
|l995h . lKettner et al dl995l). I n the phenomenological MIT 
bag model jchodos et qj|l974l ) , it is assumed that the quark 
confinement is caused by a universal pressure B on the sur- 
face of any region containing quarks. 

In the bag model, the strange matter equation of state 
(EOS) has a s imple linear form given by p = l/3(p — 4_B) 
(|witten|[l984l ). where, p is the energy density, p is the 
isotropic pressure and B is the bag constant. IDev et al 
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( 1998) developed a new model for strange stars where the 
quark interaction is described by an interquark vector po- 
tential originating from gluon exchange and a density de- 
pendent scalar potential which restore s chiral symmetr y at 
a high density The EOS formulated bv lDev et al (1 19981 ) can 
also be approximated to a line ar form p = n(p — p 3 ), where , 
n and p s are two parameters l|Gondek-Rosiriska et aj[2000l : 
IZdunik|[2OO0l ). Thus, if pulsars are actually strange stars, 
linearity in EOS seems to be a feature in the composition 
of such objects. Once the EOS is known, for a given cen- 
tral density or pressure, the Tolman-Oppenheimer-Volkoff 
(TOV) equations are then integrated to calculate the macro- 
scopic features such as mass and radius of the star. 

However, as the densities within such stars stars are nor- 
mally beyond nuclear matter density, one expects anisotropy 
to play a major role in such calculations as the TOV equa- 
tions get modified in the presence of anisotropic pressures. 
Note that among many possibilities, one of the reasons for 
the consideration of anisotropy within such stars could be 
the p r esenc e of strong electric field as recently suggested by 
lUsovl (|2004h . 

Since the initial paper bv lBowers fc Liand i| 19741) there 
has been extensive research in the study of anisotropic rel- 
ativistic matter in general relativity. The analysis of static 
spherically symmetric anisotropic fluid spheres is important 
in rela tivistic astrophysics. The investigations of iRudermanl 
||1972T ) indicated that nuclear matter may be anisotropic 
in high density ranges of order 10 15 gm cm -3 where nu- 
clear interactions have to be treated relativistically. Ce- 
lestrial bodies are not composed of perfect fluids only: 
the radial pressure is not equ al to the tangential pressure. 
iKippenhahn fc WeigertJ |l990|) showed that anisotropy could 
be introduced by the existence of a solid core or the pres- 
ence of type 3A superfluid. Anis otropy can ari se from differ- 
ent kinds of phase tran sitions (|Sokolovlll980l ) or pion con- 
densation (Sa wverlll972l) . Anisotropies in spherical galaxies, 
in the context of Newtonian gravity, have been studied by 
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iBinnev fc Tremaind d 19871). The effects of slow rotation in 
a star (jrlerrera fe Santosf 19951 ) or the mixture of two gases 
(|Letelierlll980l 7 such as ionized hydrogen and electrons, can 
be modelled as a relativistic anisotropic fluid. Strong mag- 
netic fields can generate an anisotropic pre ssure compon ent 
inside a compact sphere as pointed out by I Weber] (|l999l ). 

Our analysis depends on two key assumptions. Firstly, 
we choose a rational functional form for the mass function. 
The form chosen ensures that the mass function is contin- 
uous and well-behaved in the interior of the star; it yields 
a monotonically decreasing energy density in the stellar in- 
terior. This is a desirable feature for the model on physi- 
cal grounds. It is interestin g to observe that many of the 
solutions found previously (jChaisi fc Maharail 120051 . 120061 ) 
do not share this feature. It is important to note that the 
mass function chosen is similar to the p rofile that a ppears 
in the modelling of dark energy stars (|Lobol [2006) . Sec- 
ondly, we choose an equation of state w hich linearly relates 
the ra dial pressure to the energy density. iMaharai fc Chaisil 
(2006) have demonstrated that this simple equation of state 
is consistent in the modelling of compact matter such as 
neutron stars and quasi-stellar objects. The equation of 
state yields values of surface redshifts and masses that cor- 
respond to realistic stellar objects such as Her X-l an d 
VelaX-1 jGokhroo fc Mehralll994l ;l Chaisi fc MaharailliobH) . 
Linear equations of state are relevant in the modelling of 
static spherically symmetric ani sotropic quark matte r dis- 
tributions as demonstrated by iMak fc Harkol (2002) and 
anisotropic quark st ars admitting a conformal symmetry 
<|Mak fc Harkoll2004l ). 

We attempt here to construct a new model for compact 
stars where pressure anisotropy is present and the EOS is 
linear. Our paper is organized as follows: In section 2, we 
develop a model for anisotropic stars with an in-built linear 
EOS. The physical analysis of the model will be discussed in 
section 3. We also perform some numerical calculations to 
show the role of anisotropy on the gross features of compact 
stars. In section 4, we summarize our results and point out 
avenues for further investigation. 



pressure and pt is the tangential pressure. These equations 
may be combined to yield 



(P + PrW + 2Pr + ~(Pr ~ Pt) 







(6) 



which is the conservation equation. If the mass contained 
within a radius r of the star is defined as 



m(r) 



r p(r )ar 



we obtain the equivalent system of field equations 
, 2m(r) 
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2m(r) + p r r 3 



r(r — 2m) 



dp 
dr 



^ - -(P + Pr) 



2m(r) + r 3 p r 2A 
2r(r — 2m(r)) r 



(7) 

(8) 
(9) 
(10) 



where A = p t — p r , is the measure of anisotropy in this 
model. Equation (|10|l is the modified TOV equation in the 
presence of anisotropy. 

To make the above set of equations tractable, we choose 
the mass function in a particular form 



m(r) = 



br 3 



2(1 



(11) 



where a and b are two positive constants. The choice of 
the mass function is motivated by the fact that it gives a 
monotonically decreasing energy density in the stellar inte- 
rior. Similar forms for the mass function ha ve earlier been 
considered by iMatese fc Whitman! (fl980l ) andlFinch fc Skeal 
<| 19891 ) for isotr opic fluid spher e s. lLobd ((20061 ) for dark en- 
ergy stars, and IMak fc Harkol (|2003l ) for anisotropic fluid 
spheres. 

Substituting equation (JTTJ in equation ((H), we obtain 
one of the metric functions in the form 

1 + ar 2 



l + (a-&)r 2 ' ^ 
From equation ([3]), the energy density then can be obtained 



2 ANISOTROPIC MODEL 

We assume the line element for a static spherically symmet- 
ric anisotropic star in the standard form 



e^dt 2 + e^ r) dr 2 + r 2 (de 2 + sin 2 



2 ) (1) 

where j(r) and p(r) are the two unknown metric functions. 
Assuming the energy momentum tensor for an anisotropic 
star in the most general form 

Tij = diag (-p, p r , p t , p t ) (2) 
the field equations are obtained as 



Pr = 



7 ' e -" (l-e~ 



Pt = — — 27+7 - 7 fi H 



(3) 
(4) 
(5) 



where primes denote differentiation with respect to r. In 
equations ©-(O, p is the energy density, p r is the radial 



P = 



6(3 + ar 2 
(1 + ar 2 ) 



(13) 



so that the central density is given by p c — 3b. 

To find the other metric function y(r), in general, 
either a special fu nctional form o f the anisotropic pa- 
ramet er A(r) (see IMak fc Harkol d2003t). iDev fc Gleisei 

(2004) ) or the pressure p r (r) profile (see lChaisi fc Maharai 

(2005) ) is chosen at this point. However, since linearity 
seems to be a feature in compact star EOS (as shown by 
ISharma fc Mukherieel l|200ll ); | Sharma et a] (|2002l )). we as- 
sume a linear form of the EOS to solve equation (O. Inte- 
grating equation ([9| we get 



7 = 



p r (l + ar 2 )r 

l + {a-b)r 2 "" ' 2{a-b) 



dr+- 



■ln(l + ar 2 -6r 2 )+lnB(14) 



where, ln_B is an integration constant. To evaluate the first 
term on right hand side of equation (|14[) . we assume a linear 
EOS of the form 



p r = n(p - p s ), 



(15) 
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where ^ n ^ 1 is a constant which is related to the sound 
speed dp r /dp = n, and p s is the density at the surface r = 
s, where the radial pressure vanishes. We do not impose 
any restrictions on the tangential pressure p t . We obtain 
the other metric function in the form 

„„„2 



e 7 = B(l + ar 2 ) n (l 



or ) exp 



2(a - 6) 



where 

5nab 



k 



2na 2 



3nb 2 + ab-b 2 - be 



(16) 



(17) 



2 (a - 6) 2 

In equation (|17p we have set c = np s . Finally using equa- 
tion QlOp. we get the anisotropic parameter as 

a _ r2 (P + Pr) if-,, 2\ , ii abnr 2 (5 + ar 2 ) 



>r(l- 



i + bl 



(18) 



4(1 + ar 2 — 6r 2 ) LJ " A " ' ' ' " J (1 + ar 2 _ 

which clearly shows that anisotropy vanishes at the centre 
in this model, i.e., p r — Pt at r = 0. The metric potentials 
are nonsingular at the centre. The energy density and the 
two pressures (p r and pt) are also well behaved in the stellar 
interior as will be shown later. 

In our model, a, 6, c have the dimension of lengthT 2 . For 
simplicity in numerical calculations, we make the following 
transformations: 



6iT 



where, J? is a parameter which has the dimension of a length. 
In terms of these dimensionless parameters our results may 
be summarized as follows: 
a 1 + 5y 



P = 



Pr 
Pi 

where 

5nab 



1 + (5 - b)y 

B(l + a 2 /) n [l + (a-6)y] fc exp 

6(3 + ay) 
R' 2 (l + ay) 2 
n(p - p s ) 

Pr+A 



acy 



2(a - 6) 



k = 



3nb 2 + ab-b 2 



be 



(19) 
(20) 

(21) 

(22) 
(23) 

(24) 



2(5 - 6) 2 

and y = r 2 /R 2 . We may also write the anisotropic parameter 
in the form 

A 

A = = x 



4i? 4 (l + ay - by)(l + ay) 3 
[n(n + 1)6 2 (3 + ay) 2 - (2n + l)6c(3 + ay)(l + ay) 2 
-c 2 (l + ay) 4 + b 2 {n + 1)(3 + ay)(l + ay) 
-6c(l + ay) 3 - 4nab(5 + ay) (I + ay - by)]. (25) 
The mass contained within a radius s is given by 
bs 3 /R 2 



M 



2(1 + ds 2 /R 2 



(26) 



• The interior solution should be matched with the 
Schwarzschild exterior model at the boundary r = s, which 
gives 



a(r=s) 



1 _ 2M^~ 1 



y(r=s) _ f j _ 2M ^ 



(27) 
(28) 



• The values of a, b should be so chosen that the energy 
density p and the radial pressure p r remain positive inside 
the star. 

• The causality condition is satisfied if we restrict our 
model with ^ n ^ 1. 

• The value of A should be so chosen that < in the 
interior of the star and the tangential pressure p t 0. This 
imposes a bound on the anisotropy given by 



A< 



r 2 (p + p r ) 



4(1 + (a-6)r 2 /i? 2 ) 



[p r (l + ar 2 /7? 2 ) + 6/i? 2 



(29) 



• The model has five independent parameters a, 6, c, n 
and B. Note that R is not independent as it can be expressed 
in terms of a, 6 and 2. Two of the above parameters can be 
fixed by the matching conditions given by equations (|27|l 
and (|28p . The third parameter becomes fixed if we choose 
the central/surface density. Thus we will be left with two free 
parameters governing the geometry and EOS of the star. 

By suitably choosing values of the unknown parame- 
ters, it is possible to show that our model can describe 
realistic compact stellar objects. For example, let us con- 
sider a particular case where the central density and the 
surface density are given by p c — 4.68 x 10 15 gm cm -3 
and p s = 1.17 x 10 15 gm cm -3 , respectively. If we now 
choose, a = 53.34, 6 = 54.34 and R = 43.245 km, then 
the model yields a star of mass M — 1.435 Mq and radius 
s = 7.07 km. Note that these are the values which can be 
obtained if we use on e of th e EOS for strange matter for- 
mulated by iDev et a] (|l998f ) . The behaviour of the energy 
density, two pressures and anisotropic parameter are shown 
in Figures 1-3, re spectively . It is t o be noted here that, un- 
like the slope in IDev et al |l998) (which was n = 0.463), 
the maximum slope for a physically reasonable model in 
this case is given by n = 0.297. This suggests that if we 
want to obtain the same mass and radiu s for the central 
and surface densities as in lDev et ai (|l998h (which arguably 
described the X-ray pulsar SAX J 1808.4-3658), and include 
anisotropy which vanishes at the boundary, then the EOS 
becomes much softer, i.e. the gradient dp r /dp decreases. On 
the other hand, keeping the same slope (n = 0.463) and 
surface density (p s = 1.17 x 10 15 gm cm -3 ), will give us 
a star of mass ~ 1.96 Mq. Choosing different sets of val- 
ues will obviously give different results as shown in Table 1. 
Thus anisotropy provides a tool to consider a large class of 
degenerate states in this model. 



3 PHYSICAL ANALYSIS 

In this section, following Del gatv fc Lake! (|l998h . we impose 
some restrictions on the model to make it physically accept- 
able. 



4 DISCUSSION 

We have obtained new exact solutions to the Einstein field 
equations capable of describing relativistic compact ob- 
jects with anisotropic matter distribution and admitting 
a linear EOS. The solutions obtained are nonsingular and 
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Table 1 . Central density and mass for different anisotropic stellar 
models. We have set, s = 7.07 km, p 3 = 1.17 X 10 15 gm cm -3 
and R = 43.245 km. 



b 


a 


n 


pc 


Mass 








(gm cm -3 ) 


(M @ ) 


30 


23.681 


0.271 


2.58 


1.177 


40 


36.346 


0.277 


3.44 


1.300 


50 


48.308 


0.290 


4.31 


1.398 


53.34 


54.340 


0.297 


4.68 


1.435 


60 


59.789 


0.306 


5.17 


1.479 


70 


70.921 


0.322 


6.03 


1.548 


80 


81.786 


0.337 


6.89 


1.608 


90 


92.442 


0.353 


7.75 


1.661 


100 


102.929 


0.367 


8.61 


1.707 


183 


186.164 


0.463 


15.76 


1.962 



well behaved in the stellar interior. Un like other models 
|Mak fc Harkoll2003l ; iDev fc Gleiserl 120041 ), the form of the 
anisotropic parameter is not chosen a priori in this model. 
We have shown that the consideration of anisotropy may 
provide a wide range of consequences in the geometry and 
composition of a star. For example, this may either yield 
a softer EOS (so that the gradient dp r /dp may be differ- 
ent) or a star with different mass and radius. Our analysis 
depends critically on the choice of the mass function given 
by equation so that a linear equation of state is possi- 
ble. Note that some othe r treatments, such as the result of 
iMaharai fc Chaisil ([2006) , with a linear EOS has the unre- 
alistic feature of vanishing energy density at the boundary. 
Our model does not suffer from this defect. It is also im- 
portant to observe that in the presence of anisotropy the 
material composition, in objects such as SAX J1808. 4-3658, 
need not necessarily be quark matter because the softness 
may vary. In future work it would be interesting to investi- 
gate what other forms of the mass function are consistent 
with a linear equation of state. 
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Figure 1. Energy density plotted against radial distance. 
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Figure 2. Radial pressure (solid line) and tangential pressure 
(dashed line) plotted against distance. 
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Figure 3. Anisotropic parameter plotted against radial distance. 



